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Abstract 

We present measurements of the scale (bin-size) dependence of event-wise mean 
transverse momentum (p t ) fluctuations for Au-Au collisions at ^Js^n = 200 GeV 
simulated with the Hijing-1.37 Monte Carlo. We invert those scale distributions to 
obtain autocorrelation distributions on pseudorapidity rj and azimuth (j) difference 
variables (tja^a)- The autocorrelations have a simple structure dominated by a 
same-side (on azimuth) peak with similar widths on the two difference variables. 
With jet quenching imposed, the same-side peak widths increase slightly in both 
directions, and the amplitude is substantially reduced. The same-side peaks are 
attributed to minijets, observed in this study as local velocity correlations of charged 
particles with p t < 2 GeV/c. 
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1 Introduction 



According to QCD theory a phase transition from hadronic nuclear matter 
to a color-deconfined medium (quark-gluon plasma or QGP) may occur in 
relativistic heavy ion collisions at RHIC (jl|). Measurements of nonstatistical 
fluctuations of multiplicity n and event-wise mean transverse momentum (p t ) 
have been proposed to identify and study the QGP (0; 0; 1; § @; 0,0). Non- 
statistical fluctuations could also result from minijet production (j9j) or from 
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other dynamical processes. Minijets (hadron fragments from a minimum bias 
parton distribution) should be copiously produced in the initial stage of ultra- 
relativistic heavy ion (HI) collisions at RHIC (1 1 Of ) . A colored medium formed 
in the prehadronic phase of such collisions could cause jets/minijets to lose 
energy through induced gluon radiation (11), a dissipative process conven- 
tionally called jet quenching in the case of higher-p t partons. The properties 
of the dissipative medium would then determine specific energy loss ( )12 ) and 
other aspects of minijet modification which could be studied with fluctuation 
measurements. 



In a previous paper (|9|) we reported an analysis of (p t ) fluctuations in simulated 
Au-Au collisions at 200 GeV at a fixed scale (bin size), using Hijing-1.37 
(flil ) to study jet/minijet contributions to those fluctuations. According to 
that study minijets are the dominant source of nonstatistical (p t ) fluctuations 
in heavy ion collisions as modeled by Hijing. The (p t ) fluctuation amplitude 
quantitatively measures the extent of jet production and jet quenching in that 
model. A measurement of nonstatistic al (p t) fluctuations in Au-Au collisions at 



130 GeV by the STAR collaboration (11 4 ) was also made at a single scale (the 



(r),(j)) detector acceptance). Comparison of that STAR measurement to our 
Hijing analysis has been inconclusive concerning minijets as the fluctuation 
source. 

In this Letter we extend our fixed-scale Hijing analysis to determine the scale 
dependence of (p t ) fluctuations. By inverting that scale dependence with a 



standard numerical technique (ll5l ) we obtain two-particle correlations induced 
by minijets in the form of autocorrelation distributions. This more differen- 
tial study of fluctuations and correlations provides detailed information about 
parton dissipation and resulting minijet structure in the Hijing model and 
establishes a reference for similar studies of RHIC data. 



2 Relating Fluctuations and Autocorrelations 



The scale dependence of event-wise (p t ) fluctuations on kinematic variables 
(77, 4>) is determined by the structure of underlying two-particle correlations (ll6l ). 
Variation of nonstatistical fluctuations (measured by variance differences) with 
bin size or scale can be inverted to obtain an autocorrelation distribution on 
difference variables. The autocorrelation distribution, a projection by averag- 
ing of a two-particle distribution on a momentum space [e.g., (771,772)] on to its 
difference variable(s) (e.g., t?a = fji — 7/2), is often simply and directly inter- 
pretable in terms of physical phenomena. The actual fluctuation analysis is 
carried out on two-dimensional (2D) space (77, 0) within detector acceptance 
(A77, A0) divided into M macrobins of variable size (scale) (Sr), S(p) which are 
in turn divided into mg microbins of fixed size (e v , e^) (microbins are used in 
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the numerical integration). We restrict to ID space x (x = rj, 0) to summarize 
the method and then generalize to 2D (77, 0) for the analysis. 



A ID macrobin of size Sx contains integrated pt{Sx) (a scalar sum over particles 
in the bin) and particle multiplicity n(5x). To measure fluctuations of p t rela- 



tive to particle multiplicity n we employ {p t (8x) — n(5x) p t } / \J n(Sx) as the ba- 
sic statistical quantity and define the scale-dependent per-particle conditional 
variance of p t given multiplicity n by o~ 2 t . n (5x) = {p t (Sx) — n(8x) p t } 2 / 'n(5x) , 
where pt = pt{Ax)/n(Ax) is the inclusive mean per-particle m and the overline 
represents an average over all bins of size 5x in all events (0; [3). 

We relate that variance to an equivalent autocorrelation distribution on dif- 
ference variable 2 a = %i — %2 in two-particle space (xi,x 2 ) by the following 
argument. Each autocorrelation element A k (e x ) is an average over the k th di- 
agonal of 2D microbins in a 2D macrobin on (xi, X2) and a further average over 



all macrobins in all events (1151 ) . Given autocorrelations for object and reference 
distributions we define the difference autocorrelation AA k = A k , bj — A k>re f. 
Since the content of a ID macrobin on x is also equal to a sum over the ms ID 
microbins within it, we can relate the variance to a difference autocorrelation 
with an integral equation expressed as a discrete sum over microbins 



a lf.rX m ^) = {Pt(m s e x ) - n(m s e x ) pt} 2 /n(m s e x ) 
y {Pt(ex)-n(e x )pt} a -{p t (e x )-n(e x )p t } b 
ail ra 5 n(e x ) 



k=l- ms U l<^<m s m & ~ \ k \ V^V^b 

^ k nk AAk (Pt ■ n;e x ) 
k=i-m s n y Ak,ref(n; e x ) 

where the integration kernel K ms -k = {m$ — \k\)/ms, and the overline indicates 
an average over all macrobins of size 5x in all events. To simplify summations 
we shift the binning on x a by 1/2 bin and assume symmetry about the origin 
on xa = (k — l/2)e x , with k G [1, ms]. The kernel then becomes K ms -k = (ms — 
k + 1/2) /ms- A reference number autocorrelation is defined by Ak yT ef{n, e x ) = 
n\(e x ) = n a (e x )fib(e x ) a _ b=k (a microbin average on the k th diagonal). 

The per-hadron conditional variance difference /S.a 2 t . n (5x) = a 2 ( Sx) — o 2 H 
integrates two-particle correlations over the scale interval [0, 8x] (jlo ). The in- 
clusive single-particle variance cr| t corresponds to self pairs in two-particle 
sums. If the self pairs are excluded from autocorrelation definitions then 
o~p t:n — > Aa 2 t . n in the above equations. To obtain a correlation measure which 
is independent of the microbin size we define autocorrelation densities by 
e 2 x Ap(p t : n, ke x ) = n k /n ■ AA k (p t : n; e x ) and e 2 x p ref (n, k e x ) = A k , ref (n;e x ) 
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(we absorb the 0(1) multiplicity-ratio factor into the definition of the differ- 
ence density). We then have 



A< ;B K e x ) = 2 f e x K m , k A f Pt \ "» * ^ , (2) 
k=i ^p ref (n,ke x ) 

where Ap/ y/p re f is the desired intensive, per-particle measure of two-particle 
correlations projected onto difference variable xa- 

For the 2D scale-dependence analysis we generalize 5x — > (5r], 5(f)) to obtain 
the per-particle conditional p t variance difference Aap t . n (5n, 8(f)) in the discrete 
form 



Aa Pt :n( m s ^ n <5 = 2a Pt Aa Pt .. n (m 5 e v , n 5 e ) (3) 

— ^ 2-~, e ri e fP-rn s n s ;kl > , 

k,i=i Jp re f(n;ke v , I e^J 



with 2D kernel K mgng . k i = (nig — k + l/2)/nig ■ (ng — / + l/2)/ris representing 
the 2D binning process. This Fredholm equation can be inverted to obtain 
autocorrelations in the form Ap/ y/p re f from (p t ) fluctuation scale dependence 



on (Sri, 6(f)) ((151) 



3 Hijing Simulations 



We used event generator Hijing- 1.37 to produce minimum-bias event ensem- 
bles with ~ 1M total events of two types of Monte Carlo collision: 1) quench-off 
Hijing - jet production enabled but not jet quenching and 2) default or quench- 
on Hijing - jet production and jet quenching enabled. Study of those event 
types should help to clarify the role of jet production and jet quenching in the 
correlation structure of RHIC heavy-ion collisions. For this analysis charged 
particles with pseudorapidity \r]\ < 1, transverse momentum p t £ [0.15,2] 
GeV/c and full azimuth were accepted. Minimum-bias Hijing events were clas- 
sified according to impact parameter b (fm) to estimate collision centrality. Six 
centrality classes were defined: 80-95%, 65-80%, 50-65%, 30-50%, 15-35% and 
0-15% of the total cross section. 
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Fig. 1. Left panels: Variance difference Acr 2 t:n (GeV/c) 2 on scales (Srj, Sep) for 65-80% 
central (upper panel) and 15-30% central (lower panel) Hijing events with jet 
quenching off. Right panels: Corresponding autocorrelations Ap/^p^J (GeV/c) 2 
on difference variables (?7a><M)- 

4 Fluctuation Scale Dependence 



In Fig. CD (left panels) we present variance differences Acr 2 t:ra (<5?7, 5(f)) (GeV/c) 2 
for Hijing quench-off collisions in centrality classes 65-80% and 15-30%. The 
variance differences typically increase monotonically with pseudorapidity scale 
(Srj) but are more complex on azimuth scale (5(f)). While fluctuation measure- 
ments can certainly be compared to theoretical predictions, the underlying 
physical mechanisms are often more evident in the corresponding autocorre- 
lation distributions. 



5 Extraction of Autocorrelations 



Fig. [T] (right panels) shows 2D autocorrelation distributions inferred from the 
fluctuation scale dependence in the left panels by inverting Eq. Q (15). We 



designate the variance difference distributions on the left as data D = Act 2 



pt-n 



and the autocorrelations on the right as image I = Ap(p t : n) / J p re f (n) , re- 
lated by Eq. ([3]) in the form of matrix equation D = T I + N, where N is the 
statistical 'noise' on the measured data. The image can be inferred in prin- 
ciple by the simple matrix manipulation T _1 D = I+T _1 N. That approach 
results in an image dominated by small-wavelength noise in T _1 N, since the 
inverse matrix T _1 , essentially a differentiation, acts as a 'high-pass' filter to 
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amplify the small-wavelength components of N. To control statistical noise on 
the image a standard regularization procedure is employed (1151 ). We treat the 
bins of image hypothesis I a as free parameters in a x 2 fit and add a smoothing 
term with Lagrange multiplier a to obtain \ 2 a = l|D — TI a | | 2 + a| |LI a | | 2 . The 
first term measures the deviation of the integrated image from the data, and 
the second term, with linear operator L, measures small-wavelength structure 
on the image. The optimum choice of a maximally reduces small-wavelength 
noise while minimally distorting meaningful structure on the image. Distribu- 
tions of | |D — TIq. 1 1 2 and | |LI a | | 2 on log a, as shown in Fig. [2] (left panel), are 
used to determine the optimum a value, bracketed in that case by the vertical 
dashed lines. 
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Fig. 2. Left panel: Distribution of ||D — TI Q || 2 (dots) and ||LI Q || 2 (triangles) vs 
a for the 80-95% central quench-on data in Fig. [3] (upper-right panel). Deviation 
of ||D — TI Q || 2 from the extrapolation (horizontal solid line) indicates systematic 
error due to image smoothing. The vertical dashed lines bracket the optimal choice 
of a. Right panel: Distribution of estimated smoothing error I a — I aa on (t]a, <M)> 
also for the 80-95% central autocorrelation distribution in Fig.O Note the 10 x scale 
reduction compared to Fig. [3j 



6 Results 



Fig. [3] shows autocorrelations for three centralities (80-95%, 30-50% and 0- 
15% from top to bottom) and two Hijing configurations: quench-off (left pan- 
els) and quench-on (right panels). All autocorrelations have similar features: 
a same-side (|0a| < 71 "/2) peak and an away-side (|0a| > 7r/2) ridge. The 
same-side peak can be identified with fragments from minimum-bias partons 
(minijets), since no high-p t trigger condition was imposed. A related structure 
in RHIC data with high-p 4 trigger for angular or number (not p t ) correlations 



can be found in (1171 ) . The away-side ridge can be identified with momentum 
conservation for some combination of longitudinal string fragmentation and 
fragments from semi-hard parton scattering. There is some indication of a 
saddle structure or minimum at rj^ = in the away-side ridge. 

The variation of the near-side peak structure with centrality is limited. For 
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Fig. 3. Autocorrelations Ap/^/p re f (GeV/c) 2 for Hijing quench-off (left panels) and 
quench-on (right panels) and for centralities 80-95% (top panels), 30-50% (middle 
panels) and 0-15% (bottom panels). 

quench-off collisions there is a modest increase in peak amplitude and small 
symmetric decrease in peak widths with increasing centrality. For quench-on 
collisions those trends are reversed. The amplitude falls and the peak widths 
increase, also varying symmetrically with increasing centrality. 



7 Errors 



The statistical errors for (p t ) fluctuation measurements as in Fig. [T] (left panels) 
are typically less than 0.001 (GeV/c) 2 for all scales and centralities. Autocor- 
relation errors have two components: statistical errors and image distortion 
due to smoothing. Statistical errors on the image are found by inverting sta- 
tistical errors on the fluctuation data. The structure of that error distribution 
is typically simple and can be characterized by a single number for each au- 
tocorrelation. The smoothing error is estimated by passing inferred ima ge 1 ^ 
a second time through the integration/inversion process to obtain I aa (1151 ). 
The difference between the two image distributions provides an upper limit on 
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Fig. 4. Left panel: Correlation amplitudes vs centrality measure v (mean participant 
path length) for near-side peaks and for quench-on collisions (dots) or quench-off 
collisions (triangles). Right panel: widths for near-side peaks on 7/a (dots) and 4>& 
(triangles) and for quench-on collisions (solid symbols) and quench-off collisions 
(open symbols). 

the smoothing distortion (because some of that change could include further 
reduction of statistical noise). A typical image-error distribution is shown in 
Fig. [2] (right panel). The maxima typically occur near regions of large gra- 
dient in the autocorrelation. The per-bin systematic error on autocorrelation 
distributions is typically a few percent of maximum autocorrelation values. 



8 Model Fits 



We fit the near-side structure for six centralities and two MC states with the 
model function in Eq. (j3J): an offset plus a single 2D peak 



F = B + B ie 



VA 


Tr, 


4>A 


V2 <Trj 


+ 





(4) 



The residuals from fitted peak structures are typically less than 5% of the 
near-side peak amplitude. Resulting peak amplitudes and widths are plotted 
in Fig. H] vs centrality measure z/, which estimates mean participant path length 
in number of encountered nucleons (lla ). The lines are drawn to guide the eye, 
with the constraint that quench-off and quench-on parameters should extrap- 
olate to the same value for v = 1 (p-p collisions). Variation of the peak shape 
(exponents r) for these pt correlations was explored. The exponent values 
T v, T <j> = 1.7 ± 0.1 provided the best description for all quench-on centralities. 
The peak shapes for quench-off data were closer to gaussian, but with a sig- 
nificant nonstatistical spread in the best-fit values: r^, = 2.0 ± 0.2. There 
was no systematic difference between exponents for r/A and 4>a directions. 
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9 Discussion 



Interpretation of nonstatistical (p t ) fluctuations in terms of velocity/temperature 
correlations on (77, (f>) is described in (1191) . T he present analysis provides direct 
access to the correlations suggested in (1191 ), in the form of p t autocorrelations 
obtained from inversion of (p t ) fluctuation scale dependence. Those autocor- 
relations measure covariances of (pt) fluctuations in pairs of small 77, bins 
separated by (t7a,0a)- That statement describes the second line of Eq. (CQ). 
Nonzero p t correlations indicate that particle pt is drawn event-wise from an 
effective parent distribution which deviates, as a function of (rj, <fi) and dif- 
ferently in each event, from the inclusive pt distribution. p t autocorrelations 
reveal the two-point correlation structure on (rj, <ft) of global properties of the 
effective parent p t distribution, such as event-wise local transverse velocity 
and temperature. The near-side peaks can be interpreted as minijets, seen for 
the first time as local velocity structures in the pre-hadronic medium. 

The modest centrality dependence of the same-side peak for both quench- 
off and default Hijing simulations is notable. The centrality dependence of 
Hijing (p t ) fluctuations measured in (j^l) disagreed with RHIC measurements 
reported in (fTi ). In that earlier Hijing study (p t ) fluctuations showed little cen- 
trality dependence, whereas (p t ) fluctuations in RHIC data were observed to 



be strongly centrality dependent (Il4l ). In the present Hijing study p t autocor- 
relations provide differential confirmation of trends in the previous study, and 
also confirm our previous interpretation of Hijing (p t ) fluctuations as mainly 
due to minijets. We might expect minijet production, at least in the more 
peripheral A-A collisions, to follow binary-collision scaling: an increase nearly 
proportional to participant-nucleon mean path length (estimated by centrality 
parameter v). The near-side peak amplitude in the pt autocorrelation should 
then be appoximately proportional to z/, at least for quench-off collisions, but 
the observed increase is actually much slower. The autocorrelation changes in 
response to jet quenching are closer to expectations. When quenching is in- 
voked the near-side peak amplitude is indeed significantly reduced relative to 
quench-off, and to an extent proportional to v. Peak widths on both difference 
variables increase, symmetrically and by only 10-20%. The linear dependence 
on v of jet-quenching effects in this model is expected. Although v estimates 
mean path length for initial-state participant nucleons it is also proportional 
to the average distance through the medium traversed by an outgoing parton. 



10 Summary 



In this analysis we have measured the scale dependence of (pt) fluctuations 
for simulated Au-Au collisions in several centrality classes obtained from the 
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Hijing-1.37 Monte Carlo. We have inverted those scale distributions to obtain 
p t autocorrelation distributions. We observe distinct minijet structures as near- 
side symmetric peaks which can be interpreted as velocity/temperature cor- 
relations on (r),(j>). The p t autocorrelation structures, especially the near-side 
peak and its response to the jet-quenching mechanism in Hijing, agree quali- 
tatively with perturbative-QCD expectations for semi-hard parton scattering 
followed by subsequent gluon bremsstrahlung depending on parton path length 
in the medium. The detailed centrality dependence of minijet correlations in 
Hijing seems inconsistent with expectations for binary-collision scaling. The 
near-side peak amplitude for quench-off collisions, representing hard parton 
scattering, increases by only 30% from p-p {v — 1) to central Au-Au, where 
one would expect a 4x increase for binary-collision scaling. The near-side peak 
width variations with centrality are small and symmetric on (?7a,0a)- These 
Hijing simulation results provide a valuable reference for heavy ion collision 
data from RHIC. We have introduced a novel analysis technique which re- 
veals in heavy ion collisions new transverse-momentum correlation structures 
related to parton scattering, in-medium parton energy loss and parton frag- 
mentation. 

This work was supported in part by the Office of Science of the U.S. DoE 
under grant DE-FG03-97ER41020. 
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